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Abstract 

Interaction of intense laser beam with an inversion symmetric 
third order nonUnear medium is modeled as a quartic anharmonic 
oscillator. A first order operator solution of the model Hamiltonian 
is used to study the possibilities of generation of higher order non- 
classical states. It is found that the higher order squeezed and higher 
order antibunched states can be produced by this interaction. It is 
also shown that the higher order nonclassical states may appear sep- 
arately, i.e. a higher order antibunched state is not essentially higher 
order squeezed state and vice versa. 

1 Introduction 

A nonclassical state of the electromagnetic field is one for which the Glauber- 
Sudarshan P-function [l] is not as well defined as the probability density is [2]. 
To be precise, if P-function becomes either negative or more singular than delta 
function, then we obtain a nonclassical state. The nonclassical states do not 
have any classical analogue. Commonly, standard deviation of an observable 
is considered to be the most natural measure of quantum fiuctuation [3] as- 
sociated with that observable and the reduction of quantum fiuctuation below 
the coherent state level corresponds to a nonclassical state. For example, an 
electromagnetic field is said to be electrically squeezed field if uncertainties in 
the quadrature phase observable X reduces below the coherent state level (i.e. 
(AX)^ < i) and antibunching is defined as a phenomenon in which the fiuctua- 
tions in photon number reduces below the Poisson level (i.e. (AA^)^ < (N)) [H 
[5]. Standard deviations can also be combined to form some complex measures 
of nonclassicality, which may increase with the increasing nonclassicality. As 
an example, we can note that the total noise of a quantum state which, is a 
measure of the total fiuctuations of the amplitude, increases with the increasing 
nonclassicality in the system [6]. 

Probably, antibunching and squeezing are the most popular examples of 
nonclassical states and people have shown serious interest on these states since 
1960s. But higher order extension of these nonclassical states are only intro- 
duced in late 1980s [TlllOj. Among these higher order nonclassical effects higher 
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order squeezing is studied in detail [71 [H [TTJ [12] but the higher order antibunch- 
ing [9] is not yet studied rigorously. As a result, still we do not have answers 
to certain fundamental questions Hke: Whether higher order antibunching and 
higher order squeezing appears simultaneously or not? The present work aims 
to provide answers to this question. In order to do so, in section 2, we have 
modeled the interaction of an intense laser beam with an inversion symmetric 
third order nonlinear medium as a quartic anharmonic oscillator. A first order 
operator solution of the model Hamiltonian is also used to provide time evolu- 
tion of some useful operators. In section 3 and 4 we have theoretically studied 
the possibilities of generation of the higher order squeezed and higher order 
antibunched states respectively. We have shown that the generation of higher 
order squeezed and higher order antibunched states is possible but they may or 
may not appear simultaneously. We finish with some comments and concluding 
remarks in section 5. 



2 The model: an intense laser beam interacts with a 3rd order 
nonlinear medium 

An intense electromagnetic field interacting with a dielectric medium induces a 
macroscopic polarization (P) having a general form 

7^ = xi^ + X2^^ + X3^^^ + --- (1) 

where is the electric field and Xi is the i — th order susceptibility. Corre- 
sponding electromagnetic energy density is given by 
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where is the magnetic field. Now, if we consider an inversion symmetric 
medium then even order susceptibilities (X2,X4 etc.) would vanish. Hence the 
leading contribution to the nonlinear polarization in an inversion symmetric 
medium comes through the third order susceptibility (xa)- If we neglect the 
macroscopic magnetization (if any) then the interaction energy will be propor- 
tional to the 4-th power of the electric field. Normal mode expansion of the 
field ensures that the electric field operator for x-th mode is proportional to 
{qx + at) and the free field Hamiltonian is 

Ho=Y,^x{alax + ^) (3) 

X 

where we have chosen h — 1. 

Thus the total Hamiltonian of a physical system in which a single mode of 
intense electromagnetic field having unit frequency interacts with a 3rd order 
nonlinear non-absorbing medium is 



H = (ata+ i) + f^(a1' + a) 



_ _ (4) 
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with 

X^l={a^+a) (5) 

and 

X^-^ia^-a). (6) 

The parameter A is the coupUng constant and is a function of xs- Here we 
can note that the silica crystals which are used to construct optical fibers are 
example of third order nonlinear medium. So third order nonlinear medium 
described by the Hamiltonian ^ is also important from the point of view of 
applicability. 

The above Hamiltonian ^ represents a quartic anharmonic oscillator of unit 
mass and unit frequency. The equation of motion corresponding to ^ is 

X + X + XX^ = (7) 

which can not be solved exactly. But in the interaction picture the potential Vi 
and the time evolution operator Uj corresponding to (j4]) are respectively 

Vi{t) = exp{ia^at)\{a + a'^Y exp{—ia^at) = \{a exp(— zt) + a'^ exp{it))* (8) 

and 

Ui{t) ^l-i f Vi{h)dh + {-if I Vi{h)dh f ' Vi{t2)dt2 + (9) 

Jo Jo Jo 

Now if we assume ^ ^ 

/ Viih)dh I ' Vi{t2)dh « 1 (10) 
Jo Jq 

then we can neglect higher order terms and ( [9]) reduces to 

Ui(t) = l-if Vi{ti)dti = l-i f \{aexp{-iti) + a'' exp{iti)fdti (11) 
Jo Jo 

Thus the first order expression for time evolution of annihilation operator is 



aj{t) = uj{t)a{0)Uj{t) = a - f [6ta + 6ta''a^ + 6exp(ii) sintat2a - exp{2it) sm{2t)a'<^ 

+ 6exp(ii) sinta'l' + 2exp(— it) sinta'^] . 

(12) 

The derivation of this first order expression (for a more generalized Hamiltonian) 
is shown in detail in [I3l[l4]. Here we would Hke to note that the annihilation 
operator in the Heisenberg picture gh and that in the interaction picture a/ are 
related by qh = exp(— it)a/(i) and our solution lfT2|) is vaHd only when (fTOj) is 
satisfied. Physically this condition implies that the anharmonic term present in 
([H provides only a small perturbation. This assumption is justified because in 
a third order nonlinear medium the anharmonic constant A, which is a function 
of X3, is very small (i.e A ^ 1) [14] . 
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3 Higher order squeezing 

Higher order squeezing is defined in various ways. The definition which we have 
used in this work is by Hillary [ll] . This definition is different from that of Hong 
and Mandel [7]. Present definition is also called amplitude squared squeezing. 
According to this definition of higher order squeezing, higher order quadrature 
variables are defined as 

Y, = ^{a^^ + a^) (13) 

and 

y, = -j=(at2-a^) (14) 

From the commutation relation [Yi'l2] = i{'^N + 2) it is easy to conclude that 
a state is squeezed in Yi variable if 

(AYi)2 < (2A^ + 1) (15) 

or if, 

f = {/SYif -{2N +1) <{) (16) 
A strenuous but straight forward operator algebra yields 

AYl = {Y^) - = [2|a|2 + 1 - I [4|a|2 (2|a|2 + 3) sini sin(20 - <) - 12|a|4t sin(40) 

+ 3 (2|a|4 + 4|a|2 + l) sin^ 2t - I2\a\^ (2|a|2 + 3) sint sin(26l + t) 
- 2|a|4 sin 2i sin (2(261 -2t))]] 

(17) 

where a = |a|exp(i0) is used. Here we would like to note that since we are 
working in interaction picture we have to take all the expectation values with 
respect to the initial coherent state |q;) which is defined as a | a) = a] a). By 
taking the expectation value of N{t) — a^^^^a{t) we obtain 

{N{t)) = |ap-^ [2|ap (2|a|2 + 3) sini sin(26' - <) - sin 2i sin (2(261 - 2i))] . 

(18) 

Substituting (fT7|) and l|18p in (fTBl) we obtain a closed form analytic expression 
for / as 



/ = [-4|ap(2|a|2 + 3)sinisin(t + 26l) -4|a|4tsin(46l) 

+ (2|a|4 + 4|a|2 + l)sin2(2t)] 



(19) 



From (fT9| we can observe that / oscillates between positive and negative values 
depending upon the phase of the input coherent fight 9 and the interaction time 
t. Both of these parameters can be tuned to produce higher order squeezed 
state and to increase the depth of noclassicality by increasing the negativity of 
/. If we consider ^ = f then (fT9|) reduces to 

f = -^ [4|ap(2|a|2 + 3) sin^ t + (2|a|4 + A\a\^ + 1) sm^{2t)\ (20) 
which is always negative and thus we always have higher order squeezing. 
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4 Higher order antibunching of photons 



As we have already mentioned the higher order antibunching [9] is not yet 
studied rigorously. Using the negativity of P function [l], Lee introduced the 
criterion for HO A as 



where N is the usual number operator, A^*^*^ — N{N — 1)...{N — i + 1) is the ith 
factorial moment of number operator, () denotes the quantum average, I and 
m are integers satisfying the conditions I < m < 1 and the subscript x denotes 
a particular mode. Ba An [15] choose m = 1 and reduced the criterion of Ith 
order antibunching to 

= / m\ - 1< (22) 

or, 

(a^('+i))<(a^W)(A^.). (23) 

Physically, a state which is antibunched in ^th order has to be antibunched in 
{I — l)th order. Therefore, we can further simplify l(23l) as 

(Ar('+i)) < (nW) (N.) < (Nf^) {N^f < (Nf^) {N.f < ... < {N^r' 

(24) 

and obtain the condition for I — th order antibunching as 

d(/) = (7V('+i))-(iV,)'+i<0. (25) 



This simplified criterion l|25|) coincides exactly with the physical criterion of 
HO A introduced by Pathak and Garica [16] . 

By using (fT2| and the condition for higher order antibunching (|25l) , it is easy 
to derive that for a third order non-linear medium having inversion symmetry, 
we have 

rf(l) = [2 (2|ap + 1) sm{t ~ 26) sin(i) + |ap sin {2{t ~ 29)) sin(2t)] , 

(26) 

rf(2) = ^^^^ [sin(< - 261) sin(t) + sin {2{t - 29)) sin(2i)] (27) 

and ^ 

d(3) Mlf^L sin {2[t - 29)) sin(2t). (28) 
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Antibunching of fourth or higher order can not be observed with a first order 
solution of the model Hamiltonian. In order to study the possibilities of their 
occurrence we have to use higher order operator solutions of 

Equations l|26H28p coincides exactly with our recent result [16] which was 
reported as a special case of a generaHzed Hamiltonian. Here we can observe 
that if we choose interaction time t = 20 then d = 0. Therefore, we can observe 
higher order coherence. But for 6* = or 6* = mr, (i.e. when input is real) 
d is a sum of square terms only. So d is always positive and we have higher 
order bunching of photons. For other values of phase {9) of input radiation 
field, value of d oscillates from positive to negative, so we can observe higher 
order bunching, anti-bunching or coherence in the output depending upon the 
interaction time t. 

5 Summary and concluding remarks 

From l(20|) we know that for 6* = -I we always obtain higher order squeezed state. 
But for 61 = f, ([271) and ([28]) reduces to 

d(2) = [- sin^ {t) + sin2 (2^)] (29) 

and 

d(3)^^^sm\2t). (30) 

respectively. Here we can see that d(3) is always positive for this particular 
choice of 0, and therefore, third order antibunching will not appear simulta- 
neously with the amplitude squared squeezing. On the other hand c?(2) os- 
cillates between positive (higher order bunching) and negative (higher order 
antibunching) values. Thus we can conclude that neither second and third or- 
der antibunching appears simultaneously nor the higher order squeezing appears 
simultaneously with the higher order antibunching. Alternatively, we can state 
that, in general higher order nonclassical effects may appear separately. The 
possibility of their appearance may be tuned by tuning the phase (0) of the 
input coherent state and interaction time (i). 
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